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Abstract 



Studying published data of Maxwellian averaged neutron capture cross sections, 
we found simple phenomenological rules obeyed by the cross sections as a 
function of proton and neutron number. We use these rules to make predictions 
for cross sections of neutron capture on nuclei with proton number above 83, 
where very few data are available. 
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1 Introduction 



Theoretical descriptions of nucleosynthesis in stars rely heavily on the knowledge of capture 
cross sections of slow neutrons on nuclei. The classical model of nucleosynthesis in weak 
neutron flux is based on slow neutron capture (the s process) that occurs along a path in 
the stability valley of nuclei (see for instance, Refs. [ffj4]). The s-process evolution codes 
take into account the most important processes (those with largest cross sections) along 
the stability valley. The necessary information on the neutron capture cross sections and /3 
decay life times, needed to describe qualitatively the abundances of the s-process elements, 
is rather well known from laboratory experiments [5T48] . 



The s-process model is capable to explain the observed abundance of heavy elements 
fairly well M . The difference of observation and prediction is largely attributed to another 
process that occurs in stellar enviroment with high neutron flux, typically in supernovae. 
In such circumstances the neutron capture is very likely and neutron rich nuclei far from 
the stability valley build up very quickly due to repeated capture of neutrons. The nuclei 
produced such a way are so unstable and short-lived that experimental information about 
their capture cross sections and decay life times is not generally available. 

In a recent work we proposed a unified model of nucleosynthesis of heavy elements in 
stars 10 . That approach takes into account all possible types of production and depletion 
mechanisms and solves the whole system of differential equations numerically. The result 
of such an approach is that (instead of the s-process path) the evolution of the synthesis 
proceeds along a band in the valley of stable nuclei. The width of this band - and conse- 
quently the final abundances of nuclei - depends on the neutron flux and the capture cross 
sections on individual nuclei charactherized by both their proton and neutron numbers, 
a(Z,N), which constitutes an essential input to the model calculations. Therefore, it is 
important to learn about these cross sections as much as possible. 

In this paper, we study the general features of Maxwellian averaged neutron capture 
cross sections collected in recent compilations of data [7j[8]. In section 2 we show some 
phenomenological observations. In the following section we use those to make some order 
of magnitude predictions for the cature cross sections a(Z, N) for proton numbers Z > 83, 
where only very few data are available. Section 4 contains our conclusions. 



2 Observations 

Maxwellian averaged neutron capture cross 
nuclei and made available in public data 
review has been presented recently in Ref. 



sections (MACS) have been measured for many 
depositories. A comprehensive and complete 
|8j. Studying the available data, we can make 
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several observations: (i) although cross sections of many nuclei have been measured, there 
are still many missing, or rather uncertain data, especially for nuclei with Z > 83 (see 
Fig.[T^); (ii) the cross sections vary over very large range of values (about four orders of 
magnitude); (iii) for any fixed neutron number N the the cross section is maximal for a 
corresponding value of the proton number Z max and decreases rapidly as \Z — Z max \ increases 
(see Fig.[T]D). The last point implies that in the Z—N plain for each iV their is a unique value 
Z min (iV) where the capture cross section attains its maximal value. The existence of such a 
maximum is qualitatively easily understood: for fixed N, increasing Z starting from a small 
value of Z, the capture of an additional N stabilizes the nucleus in the strong repulsive 
Coulomb field of the protons, the binding energy per nucleon increases. However, for Z 
above some value Z max (N) the nucleus developes a neutron skin and additional neutrons 
become more and more loosely bound and capturing any further neutrons becomes less 
likely. The quantitavie understanding is certainly more complex, which however, is beyond 
the scope of the present paper. 




Figure 1: (a) MACS (at 30keV) on nuclei as a function of the proton and neutron number, 
(b) Dependence on the proton number Z of the MACS on nuclei with fixed neutron number 
N = 30, 31, 60, 61, 90 and 91 (indicated by vertical lines on Fig.pl 



If we plot the Z max (N) function then a rather simple picture emerges: it appears that 
a simple, almost linear function can describe the data, especially for small N. This feature 
becomes even more salient if we devide the nuclei into four groups according to the even/odd 
number of protons and neutrons: (i) ZmaL for Z even, N even, (ii) Zmal for Z odd, N odd, 
(iii) Zmax for Z even, N odd, and (iv) Zmal for Z odd, N even, as shown in Fig. 2 

In Fig. [2] crosses mark the values of Z max where the n-capture cross section is maximal 
for a fixed value of the neutron number N as taken from Ref. [7]. The solid lines represent 
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Figure 2: The function Z max (iV) for even-even, odd-odd, even-odd and odd-even nuclei. 
The crosses are the experimental values and the solid line represents the fit to the function 
in Eq. (TXT) . The dashed line runs through the bottom of the stability valley. 



fits of simple functions to these points in the form of 

/(*;«., M = £±£, ID 

with a x , b and c being fitted parameters, and x = ee, oo, eo, or oe. We determined the 
values of these parameters in two steps. First, we minimized the function 

214 2 

X 2 (a, b, c) = (^max(^) - f(N; a, b, c)) , (2) 
i=i 

i.e. nuclei belonging to all four groups are taken into account and all points are assumed 
to have weight a, = 1. The upper limit in each group was chosen the largest value for 
which Z max can be identified. With such a choice the we find iVmix = 50, 56, 53 and 55 
maxima in the groups of even-even, odd-odd, even-odd and odd-even nuclei, respectively 
(50 + 56 + 53 + 55 = 214). This fit gives 

a = 0.060, 6 = 0.013, c = 0.666, (3) 
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with correlation index 



1 = ./ 1 _^_ = 997 (4) 

\ ES(*-W)-z)' V 116029 ' 



i.e. the coefficient of determination is almost one, i 2 = 0.994 [Z = ^ Ylf=i ^max(^i) = 
45.5). - — ' 

In the second step we minimize the functions 

N (x) 

iV max 2 

X\a x ) = [Z^Km - f(N; a x , 0.013, 0.666)) (5) 

i=i 

separately for each group (x = ee, oo, eo, oe). These fits result are 

a ee = 2.47, a oo = -0.38, a eo = 0.07, a oe = 0.42 , (6) 

with coefficient of determination above 0.99 in all cases. 

We also exhibit the line of the stablity valley in Fig.[2j as a function of N (instead of 
the usual A = Z + N) 

_ N + g s 

stab ~ l + b 8 N°° ' 1 } 

with parameters 

a s = 0.682, b s = 0.027, c s = 0.614. (8) 

We see clearly that the highest n-capture cross sections lie above the stability valley and 
the separation grows with N. 

We can also observe regularity in the Z-dependence of the cross section at fixed N (see 
Fig.fTjo). We can extrapolate this regularity as well as the Z mSbX values to the region in the 
nuclide chart where very few data available for n-capture cross sections on nuclei (nuclei 
with proton number above 83, see Fig. [pi. 



The first observation is a simple trend in the behaviour of the function a 
o{z mg , x {N)^ . Putting cr max (iV) on a double logarithmic plot as shown in Fig. 
panel), we find that the general trend is well described by a fourth-order power function 



max 

3 



N) EE 

i (left 



<W(A0 = mb . (9) 

This general trend is slightly modulated with some oscillatory behaviour, with minima 
around magic numbers, as seen on Fig.[3]3, where the ratios of the measured cross sections 
to cr max (N) are shown. 
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Figure 3: (a) Largest neutron capture cross sections as a function of the neutron number, 
(b) Ratio of the measured largest cross sections to cr max given in Eq. (pi. 



The second observation is that if we normalize the cross sections cr(Z, N) for a fixed 
neutron number N with the largest cross section cr max (iV), then the profile of the dependence 
on the proton number is rather similar for all neutron numbers. This similarity is best seen 
if the position of the largest cross section is shifted by — Z max to zero, therefore, we define 
these normalized and shifted cross section values, 

Pn{z) = 7~j\j\ = —( \ ' ( 10 ) 

a(z max (iV)J 

for all values of N, where data are available. Then we define the average by 

N~ 



/^) = ^rZ>^)> ( u ) 

with squared standard deviation 



N=l 



N Z (N Z 



1 Nz 12 

r — [Pn(z)-p(z)\ , (12) 



N=l 



where N z is the number of available data for fixed z. This average is shown in Fig.|4} As 
seen from Fig.|4)D this function is well approximated with an almost exponential function 
in both positive and negative directions, but with different exponents. More precisely, we 
fit the logarithm of the average with quadratic functions of the form a^z 2 + hiZ + q with 
subscript of the coefficients refering to three regions in z: (i) i — 1 for z < —26, (ii) i = 2 
for —26 < z < 0, and (iii) % = 3 for < z. For i = 2 and 3 we fix q = 0. This form ensures 
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Table 1: Result of the fit to the average function p(z). 
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di 
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Ci 


X 2 /d.o.f 


1 
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1.135 


17.95 


5.29/5 
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-0.0025 


0.2658 





6.15/9 


3 


-0.0058 


-0.3948 





7.14/4 



the constraint p(0) = 1. We also require the continuity of the fitted function at z = —26. 
We measure the goodness of the fit by the weighted sum of squares 



x 2 - 



E 



[hipO 



(diZ 2 + biZ + Ci 



(13) 



summed over values of z in the three regions separately. The result of these fits is presented 
in Table [l] and shown in Fig.|4} 




Figure 4: Average of the normalized neutron capture cross sections as a function of z = 
Z — Z max . The errorbars represent the standard deviation a(z). 

Each function Pn{z) differs from the average in two ways: (i) tipically the larger N the 
wider Pn(z) (as seen on Fig.JT^i), (ii) in addition there are seemingly random fluctuations. 
The origin of the latter could be either a small physical effect, or simply error of the 
measurement: there are published values for cross sections a(Z, N) that differ by a factor 
of two. While it is difficult to consider the effect of the latter, the first effect can be taken 
into account by a simple appropiate scaling of the width of the average to those of the 
functions Pn{z), which we discuss in the next section. 
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3 Predictions 



The phenomenological observations made in the previous section can be used to make 
predictions for the order of magnitude of neutron capture cross sections in regions of the 
nuclide chart where experimental data are not available. We make these predictions in two 
steps. First we validate our procedure by comparing our predictions to measured cross 
sections. Then we use our procedure to make predictions. 

3.1 Procedure 

Our procedure relies on three pieces of information concluded from the analysis of the shape 
of ridge of Maxwellian averaged neutron capture cross sections: 

1. position of Z max as a function of the neutron number (location of the ridge top on 
the nuclide chart) obeys the simple function Eq. ([I]); 

2. values of cx max (iV) (height of the ridge for given value of Z ma , x (N)) obey the simple 
function Eq. 

3. characteristic behaviour of the average function p(z) (slope of the ridge) is as given 
by Fig.g} 

In order to predict the cross section values for fixed neutron number, we proceed along the 
following steps: 

1. Given N, find the position of Z max from Eq. (JI|, which gives two maxima, one for 
even proton numbers (ZmL) and one for odd proton numbers (Zmax). 

2. Given Z max (either Z^lx, or Zmax), position the maximum location of the average 
function p(z) to Z max . 

3. Scale the height and width of the function p(z) to the available measured data by 
performing a two-parameter fit: (i) the scale factor of the height, (ii) the scale factor 
of the width. 

The third step is hampered by the discrepancies in the measured cross section values, which 
can sometimes be quite significant as shown in Table [2] for heavy elements. Discrepancies 
exist among data for lighter elements, but generally within a factor of two [8] . 
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Table 2: Ratios of largest and smallest measured neutron captured cross sections for 
elements beyond bismuth (8). 
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3.2 Validation 

We can compare the values of the predicted cross sections to those measured experimentally 
over the regions of the nuclide chart where data are abundantly available (Z < 82). In Fig.|5] 
we show again the cross sections of Fig.flJ) together with the predicted values following from 
our procedure described in the previous subsection. Considering the simple nature of our 
procedure, the agreement between data and predictions is striking for all neutron numbers. 
Of course, the predictions rarely coincide exactly with the measurements, but the order 
of magnitude is usually correct, especially where the cross sections are large, which is the 
most important region for nucleosynthesis. Similar agreement can be observed over the 
large region of the nuclide chart where data are available. 
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Figure 5: Dependence on the proton number Z of MACS (at 30 keV) on nuclei with fixed 
neutron number N = 30, 31, 60, 61, 90 and 91: comparison of the predictions of the 
phenomenological model to measure data. 
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3.3 Predictions of unkown cross sections 

Our procedure can be used to make prediction for cross sections in regions of the nuclide 
chart where some experimental information are available, such as Z > 83. In this region 
the general trend can be fitted to the measured data to complete the ridge. With such a 
procedure we obtain cross section values shown in Table [3] We can now use those preditions 
to complete the picture exhibited on Fig. [I] The result of such completion is shown in Fig. [6] 

n-capture cross sections 




Figure 6: Ridge of MACS (at 30 keV) on nuclei as a function of the proton and neutron 
number. 



4 Conclusions 

We studied the dependence of the published MACS data on the proton and neutron number. 
We found a simple characteristic behaviour that we call the shape of the ridge of MACS in 
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Table 3: Predictions for neutron capture cross sections (in mbarns) as a function of the 
proton and neutron number for elements beyond bismuth. 
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the nuclide chart. This shape can be described by the position and height of the ridge and 
the decrease of the slope. Quantifying these characteristics, we made predictions for cross 
sections in regions of the nuclide chart where only few data are available. Such predictions 
are vital for computer programs aimed at simulating the formation of heavy elements in 
stars. 
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